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Abstract 

We examine models in which jets are responsible for the formation of 
the emission-line spiral structure in IRAS 04210+0400. The kiloparsec- 
scale radio lobes in this active galaxy appear to be related to its extended 
emission-line spiral structure. The radio structure consists mainly of ex- 
tended symmetrically bent, FRI-type lobes, which follow the emission- 
line spiral structure at their inner edge. In the central region of the 
galaxy a double radio source is observed with a separation of approxi- 
mately 1 arcsec between its components, which are extremely well aligned 
with the hotspot from which the southern lobe expands outwards. Hill et 
al. (1988) suggested a model for the emission- line spiral structure invoking 
compressed interstellar matter, which is dragged away from the original 
jet path by the rotating ambient medium. From consideration of the prop- 
agation speed of the jets and the transverse ram pressure exerted by the 
rotating environment, we exclude this scenario as a possible origin of the 
spiral structure. We favour a model in which the jets themselves are bent 
by the rotating interstellar medium and possibly follows the emission-line 
spiral arms. We present fits of the model to the observed optical spiral 
structure. High sensitivity radio observations will be required to decide 
on the nature of the peculiar spiral structure in IRAS 04210+0400. 
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1 Introduction 



The nature of the spiral structure of disk galaxies has been a widely discussed 
issue in extragalactic astronomy since their discovery. Several mechanisms have 
been put forward in explanation, but to date no general consensus has been 
reached about a single process responsible for the spiral structure of galaxies. 
The most successful theories employ density waves (Lin & Shu, 1964) or galaxy 
interactions to produce the spiral arms (e.g. Toomre & Toomre, 1972). Another 
possibility is stochastic self propagating star formation considered by Gerola and 
Scidcn (1978). Most probably, different mechanisms are at work in different 
galaxies, and these theories are complementary rather than competing. 
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On occasion, ejection phenomena have been suggested to be responsible for 
spiral arms in galaxies which did not fit into the conventional schemes. Van 
der Kruit et al. (1972) suggested that the anomalous radio and emission-line 
arms in NGC 4258 originated from an essentially instantaneous ejection of radio- 
emitting material from the centre of the galaxy in opposite directions. In their 
model, the observed radio arms are a transient phenomenon and represent the 
current position of these plasmons, which were ejected with a wide range of 
speeds. The trajectories of the individual parcels is determined by gravitational 
forces and ram pressure. Later on, several other similar models for the jets in 
NGC 4258 have been suggested which are consistent with the observation that 
the anomalous arms are bent in the same direction as the trailing spiral arms 
(e.g. Martin et al., 1989, and references therein). 

Wilson & Ulvestad (1982) described a steady state model in which a two- 
sided jet propagates roughly in the rotational plane of a galaxy. The jet is 
then bent by the ram pressure exerted transverse to the jet by the rotating 
interstellar medium. They applied this model to the S-shaped radio structures 
found on a sub-kpc scale in the Seyfert galaxies NGC 1068, NGC 4151, and the 
radio galaxy 3C 293. The result of this model are leading spiral structures in the 
galaxy. Other mechanisms to obtain S-shaped jet pairs are pressure gradients in 
the hot gaseous halo of elliptical galaxies (Smith & Norman 1981) and precession 
of the jets (e.g. Gower et al. 1982). 

Imaging, both ground based (Hill et al. 1988; Steffen et al. 1996a, 1996b) 
and by the Hubble Space Telescope (HST; Capetti et al. 1996) suggests that 
the spiral structure of the peculiar active galaxy IRAS 04210+0400 is due to a 
narrow ridge of emission line filaments. Its redshift is 0.0462 and was classified 
as a Seyfert 2 (Beichman et al., 1985), but the classification has been questioned 
by Hill et al. (1988). It shows an unsually extended and complex radio structure 
for this type of galaxies, with a total power of 2.4 • 10 23 WHz^ 1 at 20 cm. It 
has extended radio lobes, which are related to the optical spiral arms. Based 
on their observations, Hill et al. (1988) first suggested that the spiral structure 
could be due to the remnant of the interaction of the jets with the ambient 
medium. In this scenario the emission line filaments are produced during the 
passage of the jet through the interstellar medium. It is then carried away from 
its original location by the rotation of the galaxy forming a spiral structure. 

In the present paper we analyze the suggestion by Hill et al. (1988) that the 
spiral structure in IRAS 04210+ 0400 is due to the interaction of the jets with 
the ambient medium. In Section 2.2 we consider the model proposed by Hill et 
al. (1988). In Section 2.2 we analyze the model proposed by Wilson & Ulvestad 
(1982) for other Seyfert galaxies. Section || contains the discussion of our results 
and our conclusions are summarized in Section |[ 



2 The models 
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2.1 Observational evidence 

The main observational reasons for investigating a model which invokes the 
propagation of the jets through the ISM as the cause for the emission-line spiral 
structure in IRAS 04210+0400 (Hill et al. 1988, Holloway et al. 1996; Steffen et 
al. 1996 a,b) are, in brief: 

- The spiral structure is emission line dominated. 

- Radio lobes continue the optical spiral structure. 

- Radio hotspots and bifurcations of optical spiral arms coincide. 

- The inner radio double source is aligned (to within 1 deg) with the initial 
southern hotspot. 

- The jets themselves are as yet not detected. 

These are very strong arguments for the spiral structure and the radio ejecta 
in IRAS 0421 0+0400to be physically associated, rather than related by coinci- 
dence. 



2.2 Straight jets 

The first three of the crucial observations listed above could be explained by 
a bent jet following the spiral structure. However, the accurate alignment of 
the inner with the outer radio structure suggests that the jets are straight. 
Therefore, we first evaluate the case in which the jets remain unbent during 
their propagation through the galaxy. 

During the passage of a supersonic jet through the interstellar medium (ISM) 
the bowshock in front of it sweeps up and compresses the gas in its path (e.g. 
Taylor et al. 1992). Eventually, after it has moved away from the bowshock, the 
gas is entrained by the jet or merges back into the environment. What is left 
in a stationary ambient medium are filaments of emission-line gas delineating 
the path of the jet through its environment. However, in a rotating galaxy the 
gas can be dragged away from its original position along an arc of length s(r) 
following the local motion of the ISM. This situation is illustrated in Figure [j], 
where the head of the jet has reached a distance r\ since ejection from the centre 
of the galaxy. 

The shape of the spiral arms depends on the rotation curve v g {r) of the 
galaxy and on the advance speed Vj%(r) of the jet head in the ambient medium. 

We assume circular orbits of the interstellar medium at a velocity v g (r) . The 
length s(r) of the arc over which the ionized gas is dragged is then given by 




where r is the distance from the centre of the galaxy, and r\ is the present 
distance of the jet head, which propagates at the local speed Vjh(r) through the 



3 



interstellar medium. Here and in the following, the subscript '1' indicates the 
value of a quantity at the current position of the jet head r = r±. 

An important parameter of the spiral structure related to the advance speed 
Vjh is the pitch p of the spiral structure close to the jet head. It can be defined 

as 



tan 8 



ds(r) 
dr 



(2) 



where is the angle between the jet and the tangent on the spiral arm at 
the position of the jet head (see Figure [j]). If the orientation of the galaxy 
with respect to the observer is known or projection effects can be neglected, 
then this is a measurable quantity. As was discussed in Holloway et al. (1996), 
IRAS 04210+0400 is seen roughly face-on and the jets probably propagate close 
to the plane of the sky, we shall therefore ignore projection effects. The observed 
pitch p can therefore be used to estimate the advance speed Vjh{r{) ~ pv g (ri) 
at the position close to the hotspots (see below). Other parameters, like the 
shape of the rotation curve or a slow change of the jet velocity, have a strong 
influence only on the inner spiral structure. For the simple case of constant 
advance speed of the jet head Vjh, the arc-length s(r) over which the ionized 
gas at distance r moved around the centre of the galaxy (while the jet continued 
to propagate up to a distance n) is given by 

s(r) = ^(n-r) (3) 

Vjh 

The advance speed depends mainly on the jet speed Vj and the ratio r\ of the 
densities of the external medium p x and the jet plasma pj and therefore varies 
according to the external conditions. For a highly supersonic non-relativistic 
jet, as we assume is the case for IRAS 04210+0400, ram pressure arguments 
lead to a simple approximate expression for the advance speed Vjh of the jet in 
an ambient medium (e.g. Leahy 1991) 

To obtain the corresponding relativistic expression, rj is to be replaced by 
ry/7 2 , where 7 is the Lorentz factor (after Marti et al. 1994). Hence, a relativistic 
jet mimics a heavy non-relativistic jet and the advance speed approaches the 
bulk jet speed vj as the latter approaches the speed of light. Equation || can be 
combined with Equation [l] to obtain the arc-length s(r). 

2.2.1 Evaluation 

We now evaluate the applicability of this model to IRAS 04210+0400. The spiral 
structure in IRAS 04210+0400 can be roughly characterized by the pitch p of 
the spiral at the current position of the jet head. However, the exact position of 
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the jet head is difficult to determine from the observations. In our model of a jet 
crossing a shocked interface between the interstellar and intergalactic medium, 
ISM and IGM, respectively (Holloway et al., 1996, Steffen et al, 1996c), it is 
to be expected that the expansion speed of the jet material is increased in the 
IGM. Therefore, a reasonable choice seems to be an 'effective' position of the 
jet head shortly beyond the beginning of the hotspot region. We estimate the 
pitch of the spiral in IRAS 04210+0400 to be in the range 1.5-2.5 and adopt a 
value of p = 2. Knowledge of the exact value is not necessary for our analysis. 
It is only required that p is of order unity. 

Using Equations [j] & [|, the velocity of the jet Vji at the position of the 
jet head can be expressed in terms of the measurable quantity p and the other 
parameters as follows 



Extragalactic jets are expected to have very low densities compared to the 
ISM, though the range of values is not known. However, from the assumption 
of a straight jet in IRAS 04210+0400 (based on the alignment of the inner 
and outer radio structure) and the observed pitch p of the spiral arms, we can 
estimate a lower limit for the density ratio r\. The radius of curvature Rj of 
this bending is approximately given by Rj/hj « V~ lv j/ v g-> with hj representing 
the jet diameter. The jet speed is unknown. We therefore use Equation || to 
eliminate Vj and v g introducing the observable pitch p of the spiral arms. We 
find that 



at the position of the jet head. From the distance D = 4.5arcsec and size 
d < 0.5arcsec of the hotspot and a conservative upper limit of the misalignment 
with the inner double source (a < 2deg; Steffen et al. 1996b), we estimate that 
Rj/hj w D/ (2d sin a) > 130 and therefore find the density ratio between the 
ambient medium and the jet to be rj < 1/22. From this analysis it is clear that 
the jet can get away without being bent by the rotating interstellar medium 
only if it is denser than the environment. This conclusion is drawn from the 
observations far from the centre of the galaxy. Closer in, at higher environmental 
density, the density contrast between the jet and the ambient medium is required 
to be even higher. The speed of a heavy jet is approximately the same as its 
advance velocity in the ISM. Hence, in this model, the jet velocity can only be 
roughly Vj wp- v g \ ~ 400 km s^ 1 (using v g \ = 200 kms -1 ). 

The high degree of collimation of the jet inferred from the small size of the 
hotspots, suggests that the jet is in approximate pressure equilibrium with the 
environment. Then the temperature ratio between the jet and the environment 
is roughly the inverse of the density ratio, and we therefore infer a jet temper- 
ature of 



Vjl =p- v g i(l + y^T). 



(5) 





Tj w rj • T ISM <§; Tjsmi 



(7) 
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i.e. a very cold jet, possibly with properties similar to galactic jets. 

Such a heavy ballistic jet would probably not flare at the interface between 
the interstellar and the intergalactic media. To explain these observations a 
collision of the jets with large dense clouds could be invoked (Higgins et al. 
1996). However, the symmetry of the flaring of the two jets on either side, 
renders this possibility highly unlikely. From the previous considerations we 
conclude that the suggestion by Hill et al. (1988) that a pair of straight jets is 
the origin of the spiral structure in IRAS 04210+0400 is not a viable model. 



2.3 Bent jets 

As an alternative, we now analyze a model of jets bent by ram pressure of the 
rotating interstellar medium as outlined by Wilson & Ulvestad (1982). All as- 
sumptions and limitations discussed in their paper also applies to IRAS 04210+0400. 
Their main assumptions are: constant jet speed, adiabatic expansion of the jet, 
and jet propagation in the plane of rotation of the galaxy. As discussed in the 
previous section, we shall neglect projection effects. The alignment of the inner 
double radio source with the southern hotspot is assumed not to represent the 
initial jet direction, as it was the case in the previous section. 

Figure shows a schematic representation of the geometry involved in this 
model. The steady state path of the jet can be found by numerically integrating 
the equation 



Pg( r ) v g( r ) 
Pjo h v 2 

(y'y + x) 2 ' 



Pg( r ) Vg(r) 
p g (r) v 2 (r) 
(2r-i)/2r 



-i -1/2T 



y 



/2i(r+i)/2r 



(8) 



This is Equation (8) in Wilson & Ulvestad (1982). Here x and y are Carte- 
sian coordinates, r is the distance from the galactic nucleus, p is the density and 
v is the velocity. Subscripts g and j refer to quantities of the galaxy and the 
jet, respectively. The effective radius of the jet is h, tq denotes a small start- 
ing distance and Y is the adiabatic index. Initial conditions for the numerical 
integration are y = r$ and y 1 ^> 1 at x = 0. 

An important ingredient in Equation || is the rotation curve of the galaxy, 
given by v g (r). The rotation curve of IRAS 04210+0400 is unknown because 
of the disturbed line profiles and its face-on orientation with respect to the 
observer. We considered several representative rotation curves, which can be 
described with the following equation (after Binney & Tremaine, 1987): 



>A+(i)' 



(9) 



G 



where a is a scale length and £ a constant. As representative examples we choose 
a rotation curve with a fast transition from solid body rotation to constant 
velocity with a = 0.15 kpc and a slow transition with a = 4.5 kpcand £ = 2 in 
both cases. Choosing a rotation curve which is slightly rising or falling after the 
transition from solid body rotation does not make an appreciable difference. 

Following Wilson & Ulvestad (1992) we assume that the density of the galac- 
tic ambient medium remains constant within r < and drops with the square 
of the distance further out. Using parameter values typical for spiral galaxies, in 
Figure |] we compare some example trajectories obtained varying the density of 
the jet. The values used in the numerical calculations are given in Table 1. For 
each of the three examples, we varied the density of the jet in a range for which 
the 'final' direction of the jet (arbitrarily chosen at r « 12 kpc) lays within the 
radio lobes. 

The individual trajectories have been rotated such that they all pass through 
a suitable fitting point on the spiral structure. In the northern arm we chose the 
secondary maximum of emission at r = 6kpc(PA„ = 14deg). In Figure |J we 
compare the results with an [Oiii]-5007A image which shows the emission line 
spiral structure of IRAS 04210+0400 (Steffen et al. 1996c). In the frames |(a) 
and ||(b) the calculated spirals are symmetric, using rotation curves 'a' and 'b', 
respectively. They have both been fitted to the northern arm and clearly expose 
the asymmetry of the spiral arms in IRAS 04210+0400. In the southern arm the 
bending is initially stronger and seems to straighten out or even 'bend back' at 
r w 5.5 kpc (the point of strong spectral flaring; Holloway et al. 1996, Steffen 
et al. 1996c). However, the southern arm is brighter and smoother than the 
northern equivalent. This could be attributed to a higher galactic density in 
this region (rather than lower jet power, because of the stronger line emission). 
In Figure [|(c), we therefore show a specific fit to the southern arm (Case 'c') 
chosing the point of reference at r = 4.4 kpc and PA S = 166 deg (the southern 
tip of the third contour from the absolute maximum), together with Case 'b' for 
the northern arm. Here only difference between the calculations for the norther 
and southern arm is an increase of the galactic density by 30% in the south, 
sufficient to obtain the required bending of the jet. 

Note that none of these sets of parameters are unique. Similar results can be 
obtained varying the jet velocity or the density of the ambient medium. Note 
also that a similar result can be obtained by changing the density or rotation 
velocity of the galaxy, as well as the speed of the jet, since basically only the 
relative ram-pressure determines the magnitude of the bending. Shape of the 
rotation curve or the density variation as a function of distance do also influence 
the jet trajectory to some extent. Although we cannot determine unique values 
of the basic parameters involved in the model, we can at least say that the bent 
jet model is consistent with the observations to within reasonable values of the 
parameters. 

Within the framework of this model we can now predict where the underlying 
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jet should point in the inner region of IRAS 04210+0400. We suggest that 
the northern and the southern radio jets should be oriented at position angles 
PA„ w 45 deg and PA S w 135 deg with an estimated error of 10 deg, as given 
by the range of initial position angles listed in Table 1. Within the estimated 
error, these values are almost independent of the choice of parameters for the 
jet trajectory as long as it satisfactorily fits the spiral structure. 

3 Discussion 

We have examined two different candidate models for the formation of the 
emission-line spiral structure in IRAS 04210+0400 involving jets ejected from 
the galactic nucleus. We exclude the suggestion by Hill et al. (1988) of straight 
jets which interact with the ISM as a possible explanation for the observations. 
We find that a second model, considering jets bent by the rotating ambient 
medium is a better candidate. 

Fitting this model to the optical spiral, we predict that the northern and the 
southern radio jets should be oriented along a line inclined by about 45 deg to 
the north-south axis, contrary to the position angle (6 deg) of the line joining of 
the inner radio components. High sensitivity radio observations will be required 
to detect the jets connecting the inner and outer radio structures and to clarify 
the nature of the central double source and its alignment with the southern 
hotspot. 

We found that the initially stronger bending of the southern arm compared 
to the northern arm could be attributed to different densities in the interstel- 
lar medium. An asymmetric density distribution explains not only the initially 
stronger jet bending in the southern arm, but also the apparent difference in 
scale size in the two arms. In particular, the optical and radio features appear to 
be somewhat larger in the northern arm, suggesting a reduced resistence of the 
ISM to the action of the jet. The brightness of the extended northern radio emis- 
sion is also lower than its southern counterpart (Steffen et al., 1996b), another 
pointer to a reduced interaction with the environment. If IRAS 04210+0400 is 
a disk galaxy, this effect could be due to a difference in the inclination angle 
between the jet and the disk, rather than an intrinsic asymmetry of the galactic 
disk (for a discussion of the classification of IRAS 04210+0400 see Holloway et 
al., 1996). However, because of the presence of a close companion galaxy, even 
an intrinsic asymmetry due to interaction cannot be excluded. 

The model of a stationary jet bent by ram-pressure provides no direct expla- 
nation for the two strong radio components which are located within 1 arcsec 
of the centre of the galaxy. There are broad blue and red-shifted wings in the 
spectral lines which have a spatial separation similar to the radio components, 
suggesting that they are related (Holloway et al., 1996). The longslit spec- 
troscopy and also the recent HST imaging (Capetti et al., 1996) shows that the 
inner kiloparsec-scale region is very clumpy. We therefore speculate that at least 
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the northern one of the central radio components is caused by the interaction 
of the jet with a dense cloud, possibly temporarily cutting off the power supply 
to the regions further out. This would explain the reduced brightness of the 
northern arm (in the optical and the radio regime) and the northern gap in the 
emission-line intensity. Considering the parameters used in our calculations, we 
find that this can be achieved by a dense cloud with a diameter of a few tens 
parsecs following the galactic rotation at a distance of 1 kpc from the centre (as- 
suming w g (lkpc) » 100 kms^ 1 ). The size of the cloud scales with the rotation 
velocity. 

In the present form, the bent jet model assumes that the emission line gas 
is located at the position of the jets. However, since it involves a stationary jet 
configuration, emission line gas could actually be dragged away from this path. 
Therefore, it should show a rather sharp edge on the upstream and a smoother 
or filamentary transition on the side downstream of the environmental gas flow 
around the jet. 

At this time, we cannot exclude jet precession or interaction with the com- 
panion galaxy as origins for the spiral structure in IRAS 04210+0400. For fur- 
ther modeling knowlegde about the rotation of the galaxy is vital, since the 
bent jet model predicts that the spiral arms are leading, as opposed the anoma- 
lous arms NGC 4258, which are trailing, excluding the stationary ram-pressure 
bending as a viable model for this galaxy. Martin et al. (1989) suggested that 
the jet and the tunnel created in the ISM follows the rotation of the galaxy and 
is therefore straight in the corotation zone, while it stays behind further out. 

Detection of collimated jets connecting the inner and outer radio struc- 
ture and measurement of the rotation curve of IRAS 04210+0400 would en- 
able us to distinguish between models with straight jets (e.g. interaction with 
its companion or a jet following the corotation) and non-corotating precession 
or ram-pressure bending. Since the total radio power of the galaxy is only 
2.4 • 10 23 WHz -1 at 20 cm, very high sensitivity will be required to detect the 
jets themselves. The observation of a displacement between bent radio jets and 
the emission-line spiral would indicate that the situation in reality was s ome 
combination between a bent-jet model and the model discussed in Section 2.2. 



To discuss the interaction scenario for IRAS 04210+0400 in any detail, more 
spectroscopic information on the companion and a better spectroscopic coverage 
of IRAS 04210+0400 itself are needed. Deep imaging of the environment around 
IRAS 04210+0400 would be useful to search for traces of the interaction. Future 
work should also focus on the classification of the galaxy (disk or elliptical) and 
the rotation curve of this unique object. 



4 Conclusions 

We have evaluated two different candidate models for the formation of the 
emission-line spiral structure in IRAS 04210+0400 involving jets ejected from 
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the galactic nucleus. We exclude the suggestion by Hill et al. (1988) of straight 
jets which interacted with the ISM as a possible explanation for the observa- 
tions. We find that the second model of two jets bent by the rotating ambient 
medium is a better candidate. Fitting this model to the optical spiral, we pre- 
dict that the northern and the southern radio jets should be oriented at position 
angles PA„ w 45deg and PA S w 135 deg with an estimated error of lOdeg. De- 
tection of the jets on the intermediate scale between the inner and outer radio 
structures will be necessary to discriminate between the current models. High 
sensitivity radio observations will be required to search for the jets. 
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Table 1: The jet and galaxy parameters for two different rotation curves (a,b) 
are given. For each rotation curve five different jet densities were considered. For 
rotation curve 'b' two separate sets (northern and southern arm, respectively) 
are given, which are distinguished by the density of the ambient medium. PA is 
the position angle of the initial jet direction measured in degrees. Parameters 
which are the same for all calculations are: = 2.2 kpc , v m = 200 kms -1 , 
r TO =10kpc, rj = 2 kpc, and ho — 3pc. In the table distances are measured in 
kpc, velocities in 1000 kms -1 , and densities in 10~ 23 gcm -3 . 



case 



a 



PgO PjO 

30. 1.0 

30. 1.5 

30. 2.0 

30. 2.5 

30. 3.0 



PA 



a 



0.15 
0.15 
0.15 
0.15 
0.15 



92 
68 
55 
48 
43 



b 

north 



4.5 
4.5 
4.5 
4.5 
4.5 

4.5 
4.5 
4.5 
4.5 
4.5 



30 
30 
30 
30 
30 

40 
40 
40 
40 
40 



0.06 
0.08 
0.12 
0.12 
0.14 

0.06 
0.08 
0.10 
0.12 
0.14 



58 
48 
42 
38 
35 

127 
135 
141 
144 
147 



b 

south 
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Figure 1: Schematic view of the model with straight jets and a spiral produced 
by enhanced density material being dragged away from its original position by 
the rotating interstellar medium. 
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Nucleus 



Figure 2: Schematic view of the model with jets bent by the rotating interstellar 
gas. 
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Distance [kpc] 



Figure 3: Rotation curves with a sharp and a slow rise have been considered. 
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Figure 4: The observed emission-line spiral structure (Steffen et al., 1996c) can 
be reproduced by a jet bent by the ram-pressure of the rotating interstellar 
medium in the galaxy. The curves in each panel show 5 jets with different 
densities. Frames 'a' and 'b' represent calculations for different rotation curves, 
whereas frame 'c' is for the same rotation curve as 'b', but with a higher ambient 
density in the southern spiral arm. The scale of 5 kpc marked on the frames 
corresponds to 4.5 arcsec (H Q = 75 kmK)" 1 Mpc -1 ). 



